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ABSTRACT 

We have modeled the emission of H2O rotational lines from the extreme C-rich 
star IRC+10216. Our treatment of the excitation of H2O emissions takes into 
account the excitation of H2O both through collisions, and through the pumping 
of the z/2 and vibrational states by dust emission and subsequent decay to the 
ground state. Regardless of the spatial distribution of the water molecules, the 
H2O lio — loi line at 557 GHz observed by the Submillimeter Wave Astronomy 
Satellite (SWAS) is found to be pumped primarily through the absorption of 
dust-emitted photons at 6 fim in the 1^2 band. As noted by previous authors, the 
inclusion of radiative pumping lowers the ortho-H20 abundance required to ac- 
count for the 557 GHz emission, which is found to be (0.5— 1) x 10^^ if the presence 
of H2O is a consequence of vaporization of orbiting comets or Fischer- Tropsch 
catalysis. Predictions for other submillimeter H2O lines that can be observed 
by the Herschel Space Observatory (HSO) are reported. Multitransition HSO 
observations promise to reveal the spatial distribution of the circumstellar water 
vapor, discriminating among the several hypotheses that have been proposed for 
the origin of the H2O vapor in the envelope of IRC+10216. We also show that, 
for observations with HSO, the H2O lio — loi 557 GHz line affords the greatest 
sensitivity in searching for H2O in other C-rich AGB stars. 
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Subject headings: Stars: abundances — Stars: AGB and post-AGB — stars: 
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Introduction 



Th e discovery of water vapor emission in the extrer ne C-rich AGB star IR C+10216 with 
SWAS jMelnick et al.lboOli ). its confirmation by ODIN jHasegawa et al l booeh. a nd the sub- 



seque nt detectio n of other 0-beariiiiB; mol ecules like OH (IFord et al.ll2003l ). H2CO (IFord et al. 



20041 ). and C3O (ITenenbaum et al.l 120061 ). have challenged our current understanding of the 
chemistry in envelopes around C-rich AGB stars. According to standard models, essentially 
all oxygen nuclei are predicted to be locked into CO or SiO with no reservoir to form other 
0-bearing molecules (except for low abundances of species such as HCO^ in the outer en- 
velope where photochemistry is important). Therefore, the unexpectedly high abundances 
found for H2O, OH, and H2CO, indicate that several processes not included in standard 
models for C-rich environments are driving the oxygen chemistry, but the dominant water 
production mechanism is still a source of debate. 

Four distinct mechanisms have been considered as possible sources of the observed water 
vapor in IRC+10216, each one making a specific prediction for the H2O spatial distribution in 
the envelope: (i) chemistry in the inner envelope, which would imply the presence of H2O in 
the warmest and densest regions; (ii) vaporization of i cy orbiting bodies that have survived 
from the main sequence into the C-rich AGB phase (jMelnick et al.l I2OOII : iFord fc Neufeld 



2OOII ). predicting the release of H2O at intermediate radii of R, 



int 



;i -5) X 10 



15 



cm; (izz) 



grain s urface reaction s, such as the Fischer- Tropsch catalysis on the surfaces of small metallic 
grains (jWillacyl |2004| ). which predicts H2O to attain a nearly uniform abundance at radii 



larger than Rint = (1-5 — 2) x 10^^ cm; {iv) chemistry involving photodissociation products 
in the outermost regions of the envelope: a specific mechanism relying on the radiative 
association 0+H2^H20+7 has been proposed by Agundez & Cernicharo (2006, hereafter 
AC06), and predicts H2O to be present at radii higher than R^t ~ 4 x 10^^ cm. In all cases, 
H2O is expected to have a uniform abundance from Rint up to the external region where it 
is photodissociated producing OH. 

The only water line detected so far in IRC+10216, the ortho-H20 Iiq— ^loi transition 
at 557 GHz, is the ground-state transition, with the upper level at only 27 K above the 
ground rotational level, and cannot discriminate -by itself- between the processes listed 
above. Nevertheless, the launch of the Herschel Space Observatory (HSO) will allow us 
to observe other H2O lines in the submillimeter and far- infrared domains, permitting us 
to infer the region where the H2O emission arises, and will potentially favor one of the 
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proposed formation mechanisms. The Heterodyne Instrument for the Far Infrared (HIFI) 
onboard HSO will provide very high spectral resolution observations (0.14-1.0 MHz), thus 
permitting the lines to be velocity resolved. In this paper, we model the H2O emission 
from IRC+10216 to show how the H2O spatial distribution can be inferred from HSO multi- 
transition observations. Also, we explore which H2O transition provides the most sensitive 
means of searching for water vapor around AGB stars other than IRC+10216; such a search 
would determine whether the occurence of H2O in C-rich environments is widespread. 



2. Model description 



Our model for IRC+10216 assumes a stellar mass loss rate of M= 3 x 10~^ Mq yr~^. 
The gas velocity f ield, turbulent veloci ty, and H2 density and gas temperature profiles are 
those adopted by iMelnick et al.l (I2OOII ) (see references therein), and listed in Table [H A 
distance to the source D = 170 pc is assumed. 



The adopted g as temperature p rofile is an approximate fit (e.g. lGlassgoldlll996l ) to the 
results obtained by iKwan &: Linkd (Il982l ). who calculated the thermal profile through the 
circumstellar envelope of IRC+10216 from the energetic balance between gas heating and 
cooling. These calculations included gas-dust collisional heating and both adiabatic and CO 
line cooling. HCN line cooling was ignored and, subsequently, the ISO far-infrared spectrum 
of IRC+10216 s howed strong HCN rota tional emission from both the ground and excited vi- 
brational states (jCernicharo et al.lll996l ). Nevertheless, radiative transfer calculations appear 
to indicate that HCN radiative excitation through absorption of infrared photons emitted 
by warm dust dominates over collision al excitation (j Gonzalez- Alfonso &: Cernichardll999al : 
Dinh-V-Trung fc Nguyen-O-Rieull2000l ) , and in such a case the observed HCN emission would 
not involve any further cooling. 

It is shown below that the continuum emission provides the dominant excitation source 
for all H2O lines of interest, mainly through absorption of 6 /im photons in the V2 band and 
subsequent decay, and so a model for the dust emission from IRC+10216 is required. We have 
used a radiative transfer code to compute, from the balance be tween heating and coohng at 
any ra dial position in the envelope, the dust temperature profile (IGonzalez- Alfonso fc Cernicharo 
1999bl ). Our model simulates the dust properties by making use of an empirical spectral in- 
dex (3. The inferred dust parameters are listed in Table [H and the goodness of the fit, 
shown in Fig. [H ensures that the continuum radiation density is well reproduced at least at 
r > 2 X 10^^ cm, where it varies as ~ r~^. 



The radiative transfer code used to determine the H2O excitation and emission has 



-4- 



been described in iGonzalez- Alfonso &: Cernichard (119971 . Il999bl ). The calculations assume 
spherical symmetry, and the envelope is divided into 60 spherical shells to account for the 
variations in the physical and chemical conditions with distance to the star. The dust and 
H2O are assumed to be coexistent within any shell, and the non-local radiative transfer is 
simulated through rays that cross the source. The code first computes the H2O statistical 
equilibrium populations at any radial position in the envelope with the use of an iterativ e 
procedure based on an approximate Newton- Raphson operator (ISchonberg fc Kempelll986l ). 
A maximum relative variation of populations of 10^^ was established as the convergence 
criterion; a much tighter criterion was found in some test models to change line fluxes by at 
most 0.2%. Once convergence is achieved, the resulting line profiles and fluxes are calculated. 
We report the expected intensities in units of main-beam temperature [Tmb), with a half- 
power beam width of HPBW=240" to simulate SWAS measurements, and HPBW ranging 
from 12'.'5 (1703-1910 GHz) to 41" (480-640 GHz) for HSO-HIFI observations (see TabfeEl 
below). 



^Spectroscopic data for H2O were taken from the Hitran 2004 database (IRothman et al. 

2OO5I ). We included the ground and the z/2 = 1 and t's = 1 excited vibrational states, but 



the 1^3 = 1 state was found to have a negligible effect on the emission from the ground state 
lines. The z/i = 1 state was ignored because the z/i band is much weaker than the band, 
and more excited vibrational levels are also ignorable because they are weakly connected 
to the ground vibrational state and/or the continuum emission from IRC+10216 weakens 
in the near infrared. Up to 40 rotational levels per vibrational state for both ortho-H20 
and para-H20 were taken into account, but the code allows this number to be decreased in 
external regions with low excitation, with the aim of ensuring and accelerating convergence. 

The calculations simulate both collisional excitation among rotational levels of the 
ground vibrational state, and radiative excitation in all transitions. Vibrational excitation 
through collisions is negle cted. Collisional rates among the 10 lowest rotational levels for 
Tk < 140 K are taken from lPhillips. Malluendes fc GreenI (119961 ). We extrapolate these rates 
to higher Tk by multiplying the excitation rates for Tk = 140 K by exp{Aii^(l/140 — 1/Tk)}, 
where AE{K) = Eup — Eiow is the energy difference between the upper and lower levels. The 
application of this strategy appears to be justified, since the use of the lPhillips et al.l (119961 ) 
collisional rates at Tk = 140 K, together with the quoted ex trapolation, y i eld co llisional 
rates for = 20 — 120 K that reproduce the values given by iPhillips et al.l (119961 ) within 
a factor of 1.6 . Also , some estimate of the collisional rates among levels not considered by 
Phillips et al.l (119961 ) is required to simulate the collisional excitation in the innermost re- 
gions of the envelope. With this aim, we have fou nd that if the H20-He collisiona l excit ation 
rates for the 10 lowest rotational levels given by lGreen. Maluendes. fc McLean Jl993h. are 



multiplied by a global factor of 5, most of the H2O-H2 J = 1 rates given by jPhillips et al. 
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(11996) are reproduced w ithin a factor of 2.5; similarly, most of the H2O-H2 J = rates given 
by iPhillips et al.l (119961 ) are reproduced within the same factor if the Green et al. rates are 
multiplied by a factor of 1.35, the difference in reduced mass when Ho is th e collision partner 
rather than He. After applying these corrections to the lGreen et al.l (119931 ) coefficients, some 
individual rat e s are s till strongly underestimated, by 1-2 orders of magnitude, relative to the 
Phillips et al.l (119961 ) values, specially in the ca se of collisions between para-H20 Oqo and 
H2 J = 1 (see discussion in lPhillips et al.l 1 19961 ). Nevertheless, the collisional excitation of 
a high-energy level in high-density and high-temperature regions is the result of collisional 
pumping from various other levels, a nd so indivi d ual d iscrepancies are to some extent di- 
luted. In summary, we have used the iGreen et al.l (119931 ) excitation rates, multiplied by the 
factors quoted abov e, to simulate the collisional excitation among levels not considered in 
Phillips et al.l (119961 ). The adopted extrapolation does not critically affect results because, as 
shown below, radiative excitation dominates over collisional excitation in all lines of interest. 

We have generated a set of models that differ only in the inner radius of the H2O shell, 
Rint, and in the H2O abundance, X(H20). The different values of Rmt account for different 
spatial distributions of H2O in the envelope of IRC+10216, and X(H20) is derived in such a 
way that the predict ed flux of the 0-H2O Iiq— >■ loi is ~ 10~^° W cm~^ as observed by SWAS 
( iMelnick et al.ll200ll ). The ortho-to-para H2O abundance ratio is assumed to be 3:1. In all 
models, we have adopted an external radius of the H2O shell Rout = 4 x 10^^ cm, wh ere H2O 
i s assu med to be photodissociated. This value of Rout is similar to that obtained by IWillacy 
( 120041 ) and AC06 through photodissociation models of the external shells of the IRC+10216 
circumstellar envelope. The photodissociation ra dius depends on several parameters an d 
could be signiflcantly lower than the quoted value (IFord et al.ll2003l : iNetzer fc Knapplll987l ). 
but test models indicate that line flux ratios for R^t ^ 8 x 10^^ cm vary by less than 20% if 
Rout is decreased to 5 x 10^^ cm. For high values of Rint, we give in section H73] an analytic 
expression for the expected flux of the 0-H2O lio loi hne as a function of both R^nt and 
Rout, showing that results are not very sensitive to Rout as long as Rout/ Rint ^ 1- 

We have selected three of those models, called A, B, and C, that correspond to the 
proposed formation mechanisms described in section [TJ The adopted values of Rint and the 
inferred ortho-H20 abundance, X(o-H20), are listed in Table [2J Model A simulates H2O 
originating in the innermost regions of the envelope. Model B applies to both the cometary 
hypothesis and the Fischer- Tropsch catalysis mechanism. Vaporization of icy bodies predicts 
that H2O is released at a radius of (1 — 5) x 10^^ cm, depending on the evolutionary stage 
of the star after the onset of the TP-AGB phase (IFord fc Neufeldll200ll ). The formation of 
H2O due to Fischer- Tropsch catalysis on the surfaces of iron gr ains predicts t hat the H2O 
abundance attains its maximum value at (1.5 — 2) x 10^^ cm ( IWillacyl |2004| ). Therefore, 
these two models predict similar H2O spatial distributions, although the Fischer- Tropsch 
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mechanism appears to be more restrictive. Model C assumes H2O to be generated as a 
consequence of photodissociation in the outermost regions of the envelope (AC06). 



3. Model results 
3.1. The H2O lio - loi line in IRC+10216 

Figure [2] compa res the SWAS lio — ^ loi ortho-H20 continuum-subtracted spectrum (from 



Melnick et al.ll200ll ) with the line profiles obtained for models A, B, and C, with the SWAS 
beamsize. All models yield similar line profiles because in all three cases the emission region 
is small compared to the SWAS 4' beam. The line shapes obtained with the smaller ODIN 
2' beam (not shown) are also similar because, even in model C, 70% of the emission arises 
within a spherical shell of diameter 1'5 (2.3 x 10^'' cm), which is significantly smaller than 
the ODIN beamsize. The 0-H2O abundance (relative to H2) required to account for the 
observed flux, is plotted in Fig. [3] as a function of Rint', it increases slowly for low Rmt 
because the innermost regions of the envelope contribute little to the emission at 557 GHz. 
For Rint = 2 X 10^^ cm (model B), the required ortho-H20 abundance, X(o-H2 0)= 5.4x10" 



(Table El) , i s one order of raagnitu de lower than that reported previously by iMelnick et al. 



(I2OOII ) and iHasegawa et al.l (120061 ). As first noted by AC06 for the case of high Rint, this 
result is a consequence of the pumping of the z/2 = 1 vibrational state by dust emission at 6 
lim and the subsequent decay to the ground state, which dominates the excitation of H2O 
in our models. In the case of the 0-H2O 557 GHz line, the importance of mid-IR pumping 
can be shown as follows. The radiative path that dominates the pump of the lio rotational 
level is Iqi — > 1^2 lio ~^ '^21 lio- The associated pumping rate per molecule in the Iqi level 
is given by = Biu{J) ^71^72, where Biu is the Einstein coefficient for absorption of radiation 
in the Iqi — > z^2 lio line, (J) is the continuum intensity at 6.2 /xm averaged over all angles 
and integrated over the line profile, rji is the fraction of molecules in 1/2 lio that decay to 221, 
and 772 is the fraction of molecules in 221 that decay to Iiq. Ignoring fine opacity effects, and 
estimating (J) as 

= ^;^^6.2/.m, (1) 

where FQ 2^im = 2.5 x 10~^^ W cm~^ /im""^ is the observed 6.2 fim continuum flux, and r is 
the distance from the star, one obtains 

0.02 



F, ^^s-\ (2) 

where ris = r/(10^^ cm). The ro- vibrational lines are indeed optically thin in all our models. 
The 

Iqi ^ ^2 lio and Iqi —>■ 1^2 2i2 have similar opacities; for Rint ^ 2 x 10^^ cm, essentially 



- 7- 



all 0-H2O molecules are in the ground Iqi level (Fig. H^) and the radial opacities of the 
quoted lines are given by 

t:, ^ 3 X 10^2 X X ' ^ ' 



a„ / V 5 X 10-7 



17 . 

Rint J \ R'out 

which gives ~ 0.25 for model C. The fraction of 0-H2O molecules in excited energy levels 
is 20-30% at r ~ 8 X 10^^ cm and increases sharply with diminishing r (Fig. HIi); therefore, 
eq. ([3]) overestimates the opacities for Rint ^ 8 x 10^^ cm. Our models yield maximum values 
of Tr = 0.30 for the case of Rint = 4 x 10^^ cm. All other ro-vibrational 0-H2O transitions 
have excited lower levels and hence lower opacities. 

On the other hand, the collisional pumping rate from Iqi to lio per molecule in the Iqi 
level is Fc = n(B.2) Ciu(Tk), where Ciu is the collisional rate from the Iqi to the lio level, and 
can be written as 

0.01 f a„(300K) 
F^ ^ X ' 



15 



2.3 X 10-10 cm3 s-i 



X exp 1-120 ( — ^ - — ) I (4) 

where the exponential factor is an approximate fit to the variation of Ciu with T^, valid up 
to r 3 X 10^^ cm. Equations ([2]) and (jl]) show that even at ris = 1, where ~ 300 K, the 
radiative-to-collisional pumping rate ratio is ~ 2, and strongly increases with increasing r as 
a consequence of the diminishing T^. Figure Hb shows the actual values of F^ (dashed hne) 
and Fc (solid line) given by the code in model A. At ris = 1, F^ is a factor of ~ 2 higher than 
predicted by eq. ([2]) because dust and gas are mixed and (J) is higher than estimated above. 
The contribution to the excitation of other radiative paths, like Iqi 1^2 2i2 ~^ '^21 ^ lio a-iid 
2i2 ^2 loi liO) accounts for the relatively low 0-H2O abundance required to reproduce 
the observed 557 GHz flux. 

As pointed out in section [2], the inclusion of the 2/3 = 1 state has negligible effect on the 
line fluxes from the ground vibrational state, which can be shown as follows. The strongest 
line connecting the Iqi level and the z/3 = 1 state is the Iqi ^'3 2o2 transition at A = 2.63 
/im, with a pumping rate per molecule in the Iqi level given by F^s oc ^f^p ^2.6/^01; 
here g^p is the degeneracy of the upper level of the transition, A^i ^ 33 s"^ is the Einstein 
coefficient for spontaneous emission, and F2Sfim ~ 0.6 x 10"^'^ W cm~^ fim~^ is the observed 
2.6 /im continuum flux. A similar expression applies to the pumping rate in the Iqi —>■ 1^2^10 
transition at 6.2 /im [A^i ~ 10 s-^), so that the ratio of the pumping rates from the Iqi level 
into the 1^2 lio and z/3 2o2 levels is found to be ~ 55. 
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Figure [5] shows the predicted 557 GHz hne profiles as observed with HSO (HPBW=41"). 
Results are given for the telescope pointing toward the source center, and shifted 41" from 
the source center. In models A and B, the contribution to the total emission of the outermost 
regions (r > 4 x 10^^ cm, or 16") to the 557 GHz emission is low, and so the line shapes 
and relative intensity toward the offset position are typical of an optically thick, spatially 
unresolved line. Model C shows, on the contrary, a nearly flat line profile with relatively 
weak emission toward the source center because the source is partially resolved; in this case 
the intensity toward the offset position is the strongest of the three models. The predicted 
line flux and shape within the HSO beam is somewhat dependent on the assumed value of 
Rout, but significant loss of line flux in the wings of the beam is expected in model C as 
long as the source size, 2 x Rout, is larger than the beam size at 557 GHz, 10^^ cm (i.e. 

Rout/Rint > 1-25). 



3.2. Predictions for other lines in IRC+10216 

The models presented below, showing predictions for H2O lines other than the 557 
GHz one, are the same as those already described in the previous sections, with values for 
the water vapor abundance such that the 0-H2O 557 GHz line flux observed by SWAS is 
reproduced for all assumed Rint (Fig. E])- Figure [H] shows the energy level diagram of ortho 
and para-H20, and indicates the strongest lines (Ai^'.i = ±1, Ai^'+i = ±1) that lie within 
the HSO-HIFI allowed frequency range. With reference to the energy of the upper level, the 
lines may be classified into three groups: (i) 8 low-excitation lines, with Eup ^ 200 K; (ii) 3 
mid-excitation lines, with 200 K < Eup < 300 K; (iii) 6 high-excitation lines, with Eup > 300 
K {Jup > 4). Figure [7] plots the predicted line fluxes, as observed with the HSO-HIFI beam, 
versus the inner radius of the H2O shell, and Fig.[S]shows the expected line profiles for models 
A, B, and C. Line frequencies and HSO-HIFI beamwidths adopted in these calculations are 
listed in Table [31 

The overall excitation is, as in the case of the 557 GHz line, dominated by the pumping 
of the 1^2 state through absorption of 6 /im dust-emitted photons. Nevertheless, absorption 
of continuum photons in far-infrared pure-rotational transitions also plays an important role 
in the excitation of some mid- and high-excitation lines. For example, absorption of 45 
fj,m and 58 /im photons in the 523— >4i4 and 422— >-3i3 lines, contributes significantly to the 
excitation of the 523— >5i4 and 422— >-4i3 lines, respectively. On the other hand, we find that 
the effect of collisions is ignorable in all our models. We have run model A by quenching 
off all collisions among rotational levels, and found that the line fluxes vary by less than 
4%. CoUisional rates much higher than estimated in our extrapolation would be required 
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for collisions to compete with the radiation field. We conclude that the excitation of H2O 
is dominated by the radiation field in all lines and models and our results are dependent 
on neither the H2 density nor the profiles, but rely on the adequate simulation of the 
radiation field generated by dust. 

Line fiuxes decrease with increasing Rint because (i) owing to the decrease of the overall 
H2O excitation with increasing distance from the star (Fig. H^), the emissivity of all lines 
decreases with increasing R^nt more steeply than that of the lowest-lying lio — >■ loi transition, 
and since the SWAS fiux of the latter line is kept constant for all Rmt, all other line fiuxes 
decrease; [ii) for high Rint and low-excitation lines, the emission becomes spatially resolved, 
and hence line fiux is lost beyond the wings of the HSO-HIFI beam. The latter effect becomes 
quite pronounced for high-frequency lines, and explains why the 2i2 — > Iqi line (HPBW=14"5, 
corresponding to a radius of ~ 2 x 10^^ cm) is predicted to be in absorption against the 179.5 
fim continuum emission for Rint ^ 4 x 10^^ cm (Fig. [8]). 

Lines belonging to different groups trace different ranges of Rint- The low-excitation 
lines will allow us to check whether H2O is formed as a result of photodissociation products. 
If H2O is formed in the outermost layers of the envelope, with Rint ^ 2 x 10^^ cm, the 557 
GHz line is expected to be the strongest of all lines. Most lines would be undetectable in 
such a case, given the expected sensitivity of HSO-HIFI (Table [3]). 

The observation of the high and mid-excitation lines would check whether H2O is formed 
in the innermost layers of the envelope. The fluxes predicted for the high-excitation lines 
decrease by one order of magnitude or more for 4.5 x 10^^ cm < Rint < 2 x 10^^ cm; however, 
these high-frequency lines are hardly detectable with the expected sensitivity of HSO-HIFI. 
Also, the predicted fluxes are somewhat uncertain because they depend on the details of 
the radiation fleld in the innermost regions of the envelope, which in turn depend on the 
distribution and emissivity of dust in these regions. Nevertheless, some mid-excitation lines 
are still sensitive to low values of Rint, as, for example, the 812 -^Sos and 821 — >3i2 transitions. 
These lines, as well as other low-excitation lines, will also probe intermediate values of Rint- 
The fluxes of the 3i2^3o3 and 812^221 lines vary by a factor of ~ 2.5 from Rint = 2 x 10^^ 
cm to Rint = 4 X 10^^ cm, and the flux of the 3o3— ^2i2 line varies by a factor of ~ 2 in the 
same Rint interval. 

The expected line proflles (Fig. [8]) display a variety of shapes, from parabolic (optically- 
thick, spatially-unresolved lines) to triangular (high-excitation lines), double-peaked (spatially- 
resolved lines), top-hat (optically-thin, spatially- unresolved lines) and P-Cygni (spatially- 
resolved, ground-state lines) proflles. The greater uncertainty concerns the predicted shapes 
of the high-excitation lines, as they are formed in the acceleration region. 
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3.3. HSO observations of C-rich AGB stars other than IRC+10216 



In this section we explore the prospects for detecting any of the low-lying ground vi- 
brational state H2O lines, lio— >loi at 557 GHz, In— >Ooo at 1113 GHz, and 2i2— >loi at 
1670 GHz, in C-rich AGB stars other than IRC+10216. Figure [9] shows that, for a source 
similar to IRC+10216, the line flux ratios depend for high Ri^t on the distance to the source, 
as a consequence of the spatially extended emission and of the different HSO/HIFI beams 
available for the three lines (Table [3]). For Rmt ~ 4 x 10^^ cm, the 1670 GHz line will only 
be spatially unresolved if the source is located at more than 2 kpc. We will restrict in the 
following to models with Rint = 2 x 10^^ cm, which predict spatially unresolved emission 
in all lines at essentially any distance. This value of Rint corresponds to model B in the 
previous sections, which assumes H2O to be formed or released at intermediate radii, and 
whose results for the low-lying lines considered here are also very similar to those obtained 
for the lowest Rmt (model A, see Table |2]). 

The role of the mid-IR continuum emission in the excitation of H2O in IRC+10216 
should be also evaluated in other C-rich stars. To this end, we show in Fig. [TOb the 6.3 
mid-IR flux, corrected for the distance, versus the mass loss rate for a sam ple of C-rich AGB 



stars. Most M and D values have been taken from the compilation of iGuandalini et al. 



( I2OO6I ). Distance estimates were mo s tly in ferred from HIPPARCOS astrometric measure- 
ments by iBergeat. Knapik. fc Rutilyl (120021 ) ; the mass loss ra tes, derived frorn millimetric 



molecular observations, were extracted from th e literature (e.g..lLoup et al 



quently corrected for the updated distances by lBergeat &: Chevallierl (120051 ). Filled symbols 



19931 ) and subse- 



in Fig. [TOb indicat e sources for wh ich the 6.3 /xm flux has been directly measured from 
ISO/SWS spectra Jsioan et al.ll2003l in. The ISO/SWS spectra has also allowed us to deter- 
mine the 6.3-to-8.8 /^m flux density ratio as a function of M. This ratio, ranging from 1.25 
for the highest M to 3.5 for the lowest M, has been used to estimate the 6.3 fim fluxes, 
shown with open sym bols, for those sources not observed by ISO/SWS but for which the 8.8 
fim flux is available (IGuandalini et al.l l2006l : iMonnier et al.l 1 19981 : iLopez et al.l 119931 ) . The 
6.3 fim emission is responsible for the excitation of H2O through the pumping of the z/2 = 1 
state, and is thus a measure of the impact of the mid-IR emission on the excitation of H2O. 
The plot shows that IRC+10216 has a high mid-IR intrinsic luminosity relative to other 
C-rich stars with similar M. In order to simulate the general trend Lg.s/im — M found for 
the bulk of sources, we have generated dust models that differ in both M and the stellar 
luminosity (L^,). For the aim of simplicity, only has been varied as a function of M, while 



-'^An Atlas of fully processed ISO/SWS spectra is publicly available at 
http: / /isc. astro.cornell.edu/ ^sloan/library/ swsatlas/atlas.html 
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the rest of dust parameters remain equal to those adopted for IRC+10216 (Table [T]). The 
grey line in Fig. [TOk shows the fit found with = 3.5 x 10^ Lq for M= 1.2 x 10"'' Mq yr~\ 
and = 8 X 10^ for M= 3 x 10"^ M© yr"^ 

Figure [TOb shows the 6.3 /xm flux density, corrected for the distance and divided by M, 
as a function of M. At any radial position of a given source, the mid-IR radiation density is 
expected to be proportional to F^ ^^^D"^ and the H2 density proportional to M, so that the 
quantity in the ordinates of Fig. [TUb is proportional to the radiative-to-collisional pumping 
rate ratio. If the kinetic temperature profile is not very different from that adopted for 
IRC+10216, it is expected that H2O in those sources with Fg.s^m-D^/M similar to or higher 
than the value associated with IRC+10216 will be radiatively pumped. Figure [TOb shows 
that radiative pumping is expected to dominate over coUisional pumping at least for the 
bulk of sources with M < 10~^ M© yr~^. 

We have used the above dust models to compute the expected H2O emission as a function 
of M. The key assumption for the H2O models displayed below is that the H2O outflow rate 
(i.e. the total number of H2O molecules in the envelope) is independent of M and equal 
to the value found in IRC+10216 (which was required to explain the 557 GHz emission). 
That assumption is appropriate if the origin of the water vapor is the vaporization of icy 
bodies, because the amount of comets orbiting a star is expected to be independent of M, 
whereas any other formation mechanism involves chemical reactions whose volume rates are 
proportional to the square of the density. In the cometary hypothesis, the release of H2O to 
the outflow is still a function of the stellar luminosity, which is expected to increase with M, 
but also depends on other p arameters such as the evolutionary stage of the star after the 



onset of the TP-AGB phase ( jFord fc Neufeldll200ll ). and thus a correlation between the H2O 



outflow rate and M is uncertain but probably less pronounced than for any other formation 
mechanism. 

Another important model assumption is the parametrization of the photodissociation 



radius for H2O molecules. Rout-, as a function of the mass loss rate. FoUowing lNetzer &: Knapp 



(119871 ). we have assumed that Rout scales as M° '', but normalized the values of Rout in such 
a way that Rout = 4 x 10^^ cm for M= 3 x 10"^ M© yr"^ (Willacy 2004, AC06). As a result 
of this different normalization, our valu es for Rprnt are one orde r of magnitude larger than 



the radius of the OH shell calculated by Netzer fc KnappI ( 1987 ). For high M, our resulting 



values for Rout are much larger than the assumed Rmt = 2.1 x 10^^ cm, and results are not 
sensitive to the adopted value for Rout- However, for M< 10~^ Mq yr~\ the finite size for 



the H2O shell (e.g. Rput/R int S, 10) starts to 



ower the expected line fluxes, in particular that 



of the 557 GHz line. Since iNetzer &: KnappI (119871 ) obtain values for Rout that are still much 



lower than ours, one should consider with caution the H2O fluxes predicted for M< 10" 
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Mq yr~^. Nevertheless, we note that the photodissociation rate in the unshielded interstellar 
radiation field (5.9 x 10~^° s~^) yields, for a 15 km s~^ fiow, a lengthscale of ~ 2 x 10^^ cm. 
For Rint = 2 X 10^^ cm, this implies a minimum value for Rout of ~ 4 x 10^^ cm for any M 
however small, unless the radiation field is stronger than the average value. For the lowest 
M= 1.2 X 10"'^ Mq yr"^ considered here, the adopted Rout is 8.4 x 10^^ cm. 

Figure [TT] shows the expected fiuxes of the three ground-state lines of H2O as a function 
of a) M, and b) the 6.3 /im continuum fiux; results are given for a source at D = 0.5 kpc. We 
find that the excitation is dominated by absorption of mid-IR photons; only for the highest 
M does coUisional excitation account for at most 20% of the predicted fiux in the ground- 
state lines. Therefore, a tight correlation between line fiuxes and the 6.3 /im continuum fiux 
density is found in Fig. [TTb . where fits of the form F oc F^^^^ are displayed. The slope 
6 is ~ 1 for the 1670 GHz line, but somewhat lower for the 557 and 1113 GHz lines (0.87 
and 0.94, respectively) due to slight opacity effects in the vibrational lines responsible for 
the radiative pumping. The slopes would approach to unity for a microturbulent velocity 
higher than the assumed value of 0.65 kms~^. For low M, the data points are slightly below 
the fitted lines as a result of the H2O photodissociation effects mentioned above; stronger 
departures from the fitted lines are expected for lower values of Rout- 

According to the current estimates of the HSO-HIFI sensitivities (Table [3]), the 557 GHz 
line will be detected in 1.5 hours at more than 5(T level if the line fiux is higher than 1.3 x 10"^^ 
W cm~^ (after coadding the two polarization modes). From Fig. [TTb . this sensitivity limit 
translates into a 6.3 /xm fiux density of 1.7 x 10~^^ W cm~^ /xm"^ aX D = 0.5 kpc, which can 
be in turn adapted to any other distance provided that both the fine and continuum fiuxes 
vary as D~^. Figure [TTb shows the 6.3 /xm fiuxes for the sample of C-rich AGB stars, with the 
solid line indicating the 6.3 /im fiux density required to detect the 557 GHz line according to 
the quoted sensitivity limit. All sources above the solid line would he detected with 1.5 hours 
of observing time in the 557 GHz line if the H2 outflow rate were equal or higher than the 
value found in IRC+10216, and if the photodissociation radius of the H2O molecules remains 
sufficiently large relative to Rmt- The sensitivity limits for the 1113 and 1670 GHz lines are 
also displayed, and show that the detection of these lines require significantly higher 6.3 /im 
continuum fiuxes. We conclude that the 557 GHz line provides the most sensitive means of 
searching for H2O in C-rich AGB stars other than IRC-l-10216 with HSO. 
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Discussion 



4.1. Model uncertainties 



Our models for lRC+10216 predict that, whatever the region where H2O is formed or 
released to the outflow, H2O is radiatively excited, which implies that model results are 
independent of the gas temperature and density p rofiles. This re s ult re lies on the extrapo- 
lation to higher Tk of the coUisional rates given by iPhillips et al.l (119961 ) , which will require 
further confirmation from new estimates of H2O-H2 collisional rates at higher temperatures. 
Nevertheless, we find that collisional rates would have to be about one order of magnitude 
higher than estimated to compete efficiently with the radiative rates in the innermost regions 
of the envelope, which seems somewhat implausible. 

A relatively uncertain parameter in our models is the assumed distance to IRC+10216, 
D = 170 PC. It has been proposed that the dis t ance may be substantial ly smaller, D = 
100 — 150 pc (IZuckerman. Dyck. fc Claussenlll986l : iGroenewegen et al.lll992l ). At 170 pc, the 
inferred stellar luminosity is L^, fa 2 x 10^ L©, which results in L^, fa 10^ Lq at D = 120 pc, 
a value more similar to the inferred l uminosities of other C-rich AGB stars with similar M 
(ISchoier. Olofsson. &: Lundgreru |2006| ) . U D = 120 pc, both M and the radiation density at 



6 fim are a factor of 2 lower than assumed, so that the radiative-to-collisional pumping rate 
ratio still remains unchanged. Since the closer proximity and weaker 6 fim radiation density 
have opposite effects on the H2O outflow rate required to match the SWAS 557 GHz flux, 
the latter will decrease by less than a factor 2, and thus the expected H2O abundance will 
be a factor of < 2 higher than in Tabled Therefore, we conservatively estimate an 0-H2O 
abundance in the range (0.5 — 1) x 10~^ for Hint = 2 x 10^^ cm. Concerning the line ratios 
to be observed by HSO, we expect values similar to those obtained at D = 170 pc if the line 
emission remains unresolved, i.e. for low values of Rmt- For high values of Rint, beam effects 
will be more important a.t D = 120 pc, diminishing the low-excitation line fluxes relative to 
the 557 GHz line. 

A potentially more important source of uncertainty is the assumption of spherical sym- 
metry in our models. Both the molecular and dust emission from IRC+10216 show evidence 
for departures from a smooth distribution in a spherically symmetric outflow; incomplete, dis- 
crete shells or arcs and clumpy structures are instead observed on a wide range of distances to 



the star (e.g. iLucas fc Guelinlll999l : iMauron &: Hugginslll999l : iFong. Meixner. &: Shahl 12003 



Leao et al.ll2006h . Even more important, the OH line shapes observed in IRC+102 16 strongly 



sugge st an asymmetric distribution of OH in the outer regions of the envelope (IFord et al. 



20031 ). thus also suggesting an asymmetric distribution of the parent H2O molecule. These 
asymmetries may alter to some extent the H2O hue flux ratios calculated with the use of 
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our spherically symmetric approach. While in spherical symmetry the emission from any 
line is isotropic, in filamentary structures or slabs of velocity-coherent gas the optically 
thick lines radiate preferentially in th e direction perpendicular to the two faces of the sheet 



( lElitzur. HoUenbach. &: McKeelll989l ). whereas the emission from thinner lines will approach 
a more isotropic behavior. Future HSO observations will show the importance of the ob- 
served morphological complexity on the line fluxes by showing whether the different H2O 
line flux ratios are consistent with a single value of Rmt, or indicate a range of values. 

Finally, the models assume a water shell with uniform H2O abundance and sharp inner 
and outer radii; however, a finite abundance gradient obviously takes place at both the H2O 
formation and dissociation regions, and variations of the H2O abundance across the shell 
are also possible. These effects may also alter to some extent the expected line flux ratios. 
Uncertainties in the outer radius of the H2O shell may also affect the expected fluxes from 
stars with low mass loss rates. 



4.2. Water formation at inner or intermediate radii 

Both the cometary and Fischer- Tropsch catalysis hypothesis predict H2O to be released 
or formed at intermediate radii of a few x 10^^ cm, and although no specific mechanism has 
been proposed for H2O formation at the innermost regions {Rint < 10^^ cm), this possibility 
cannot be rejected. The models shown in section [3l2] will permit us to discriminate, from HSO 
observations of IRC-l-10216, if there are significant amounts of H2O in the innermost regions 
through the observation of mid-excitation H2O transitions. More difficult will be a priori to 
discriminate between the cometary and Fischer- Tropsch catalysis propositions. Both predict 
similar values for Rinf, nevertheless, if Rmt were found to be significantly higher than 2 x 10^^ 
cm, the release of H2O from comets could be favoured, unless some additi onal mechanism 



were found to shift Rmt, within the Fischer- Tropsch catalysis framework (jWillacyi l2004j ) . 
outwards in the envelope. The search for water emission at 557 GHz in C-rich AGB stars 
other than IRC+10216 will also favour one of the two hypotheses: since IRC+10216 is at the 
high end of mass loss rates from C-rich stars, and the efficiency of Fischer- Tropsch catalysis 
to form H2O molecules is expected to decrease with diminishing M, one would expect in 
such rate of detection significantly lower than that quantified in section 13.31 for the 

cometary hypothesis, and one would expect a pronounced decline of the H2O abundance 
with diminishing M. 



The ISO/LWS spectrum of IRC+10216 (jCernicharo et al.lll996l ) shows an emission fea- 



ture at 179.5 fim, coincident with the wavelength of the 2i2^1oi 0-H2O transition, with a 
flux of ~ 8 X 10~^° W cm~^. This flux is very similar to that computed for the quoted line 
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in model B {Rint = 2.1 x 10 cm). However, the far-infrared spectrum of IRC+10216 shows 
vibrationally excited rotational emission of HCN, and the combined emission of the vi = 1 
J = 19 ^ 18 and 1/3 = 1 J = 19 — > 18 HCN lines, both emit ting at 179.5 /xm, is exp ected to 
be also comparable to the measured line flux a t 179.5 ^m (ICernicharo et al.lll996l ). Given 
the uncertainties inherent to the HCN model in Cernicharo et al. ( 19961 ). where the excited 
vibrational states are assumed to be thermalized, and given the high density of spectral lines 
in the IRC+10216 far-infrared spectrum, which makes it difficult to establish the contribu- 
tion from the Vi = 1 and z/3 = 1 rotational lines to the spectrum, the relative contribution 
from HCN and H2O to the observed spectral feature is quite uncertain. On the other hand, 
the expected flux of the 3o3^2i2 line at 174.6 /xm in model i? is ~ 4 x 10~^° W cm~^, 
below the 3-0" upper limit derived for that line from the ISO/LWS spectrum. These con- 
siderations suggest weakly that H2O in IRC+10216 is formed or released at radial distances 
Ri.t > 2 X 10 



15 



cm. 



Our models show that the required H2O abundance for Rint in the range (2 — 5) x 10^^ 
cm is (0.5 — 1) X 10"'' relative to H2, which makes the requirements previously demanded for 
the comet ary and Fischer- Tropsch ca talysis hypothesis to work in IRC+10216 less restrictive 
kord k NeufeldlboOll : Iwillacvl [200I . The water outflow rate is (0.5 - 1) x 10"^ M® yr'^, 
which yields, within the framework of vaporization of icy bodies, a require d total initial ice 
mass of (0.5 - 10) M® for M(H20)/Mo(ice) in the range 10"^ - 10~^ yr"^ flPord fc Neiifeld 



2OOII ). On the other hand, the Fischer- Tropsch catalysis on metallic grains will require a 
density of iron grains relative to total H nuclei of n ^fFel/rin = (0.5 — 1) x 10^^'^ to explain 
the observed H2O 557 GHz emission (jWillacyl l2004j ) . 



4.3. Water formation in the outermost layers 

Multitransition HSO observations of H2O in IRC+10216 will easily establish whether 
or not water is formed in the external layers of the envelope; the predicted fluxes in Fig. [7] 
and the HSO-HIFI sensitivities in Table [3] indicate that only the lio^loi 0-H2O line, and 
possibly the In — > Oqo and 2o2 In P-H2O lines, are detectable in one hour of observing time 
for Rint > 4 X 10^6 cm. 

For high Rint-, the expected flux in the 557 GHz line can be obtained analytically. Since 
essentially all 0-H2O molecules are in the ground Iqi level (Fig. H^), the 0-H2O Iiq-^Iqi line 
flux is derived from the radiative pumping rate given in eq. ([2]) after multiplying by 1.35 
to account for the two radiative pumping routes that are relevant at high distances from the 
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star (Fig. Hb): 

,'4 X 10^"^ cm 

F(lio ^ loi) = 1.4 X 10-^^ X 



■out 



^(Rout \ / X(o-H,0) \ W 
J V 5x10-^ J 

where X is the abundance relative to H2. Equation ([5]) is independent of the assumed distance 
D to the star because, in order to match the observed mid-IR continuum, oc (eq. [1]); 
it assumes a water shell with sharp edges at Rmt and Rout and constant H2O abundance; it 
ignores slight opacity effects in the ro- vibrational lines (section 13. ip as well as beam effects 
(both slightly raise the required abundance), and overestimates by less than 20% the 557 
GHz line flux obtained in mo del C. In the limit Rout — *■ 00, and using -F(lio —>■ loi) = 10~^° 



W cm-2 flMelnick et al.ll200lh . eq. ([5]) gives 



X(H.0) > 2.4 X 10- X (^^^) . (6) 

where x(H20) is here the H2O (ortho+para) abundance relative to H nuclei. Equation [6] 
gives the sharp-inner edge, lower limit for the H2O abundance required to account for the 
observed 557 GHz line flux, as a function of Rmt- AC06 have reported an abundance of 
x(H20) ~ 10^'' to account for the observed 557 GHz line flux in IRC+10216. However, 
the quoted abundance would imply an inner radius of Rmt ~ 2 x 10^^ cm, but the H2O 
abundance shown by AC06 (their Fig. 7) decreases sharply at r < 4 x 10^^ cm. Based on 
detailed modelling, we indicate that the x(H20) proflle given by AC06 has to be shifted 
up by a factor of ~ 2.5 to account more accurately for the 557 GHz line flux measured by 
SWAS. 

In the model proposed by AC06, atomic oxygen is produced in a shell by the photodisso- 
ciation of CO (and particularly the ^^CO isotopologue, which shields itself far less effectively 
than ^^CO.) Since the temperature is low (~ 10 K) within the shell where the atomic oxygen 
abundance is signiflcant, the neutral-neutral reaction sequence 

O + H2 OH + H 

OH + H2 ^ H2O + H 

is very slow -and therefore is negligible as a source of H2O- the flrst reaction being en- 
dothermic and the second -although exothermic- possessing a substantial activation energy 
barrier. 

AC06 therefore proposed the radiative association reaction 

O + H2 ^ H2O + 7 
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as an alternative source of H2O. To match the SWAS- and ODIN-observed water hne fluxes, 
AC06 had to posit that this reaction is relatively rapid at low temperature, with a rate 
coefficient kra ~ 10~^^cm^s~^. According to our estimation above, kra has to be a factor 
^ 2.5 higher than this value to account for the measured Iiq— >loi flux. Nevertheless, 
and even with the use of the kra estimation given by AC06, we find that this large reaction 
rate coefficient for the radiative association of O and H2 is inconsistent with observations 
of H2O and OH in at least one translucent molecular cloud, for which sufficient data exist 
to disipate any signific ant ambiguity: obser vations of the cloud along the sight-line to HD 
154368 -carried out by ISpaans et al.l (119981 ) with the use of the Goddard High Resolution 
Spectrograph (GHRS) on the Hubble Space Telescope (HST) - yield a 3 a upper limit on 
the water abundance that lies almost two orders of magnitude below the value that would 
obtain were kra as large as 10 
km ~ 2.5 X 10"^^ cm^ s~^ 



cm'^ s ^. The factor of discrepancy would rise above 150 for 



The best-fit model for the HD 154368 cloud obtained by ISpaans et al.l (119981 ) posits a 
plane parallel cloud of total visual extinction Ay = 2.65 mag in which the density of H 
nuclei is nn = 325 cm~^ (this is probably a lower limit at the cloud ce nter, as the density 



inferred from the CO J 



0/3 



2 ratio is 



10^ cm ^; 



see 



van Dishoeck et al. 



19911) a r id the external ultraviolet radiation field is 3 times mean interstellar value given by 
Draind (119781 ): luv = 3. Under these conditions, the destruction of H2O is dominated by 
photodissociation at a rate C = 5.9 x 10^^°s~^ x /[/y x (exp{ — 1.774yi} + exp{ — 1.7Av2})/2 
jLe Teuff. Millar, fc Markw"i3l200ol lN. where Ayi is the visual extinction to one cloud surface 
and Av2 = 2.65 — Ayi is the extinction to the other. At the cloud center, the water 
photodissociation rate is 1.9 x 10^ 



-10 s-i. 



For a radiative association rate of kra, the ratio of water vapor to atomic oxygen is 



therefore given by the expression 



n(H20)/n(0) = KanufHi/C, 



where /hj = ^(H2)/nH. At the cloud center, the cloud is almost fully molecular, with 
~ 0.5, and the above equation yields n(H20)/n(0) = 8.7 x 10~^ (kra/lO~ ^^ cm^ s~M. In 



Fig. [T2|, we show the predicted ?7,(H2)/r;,H ratio for the best-fit ISpaans et al.l (119981 ) model, 
together with the n(H20)/'n,(0) ratio that would result if kra were 10~^^ cm^ s~^. The model 



^We note that the factor of 1.7 in the exponential, widely used in the literature, yields H2O photodis- 
sociation rates at th e midplane of the plane-parallel cloud that are higher than the values reported by 



Roberge et al.l (|l99l[ ) by factors of 2.2 and 10 for Ay* = 1 mag and Ay* = 10 mag, respectively; therefore. 



our estimation for the -/V(H20)/A^(0) ratio predicted by the radiative association of O and H2 is probably 
a conservative lower limit. 
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predicts hydrogen to be predominantly in molecular form, in agreement with results by 
Snow et al.l (Il996l ). Averaging the ?T,(H20)/r;,(0) ratio over the entire sight-line, we obtain a 
column density ratio N{}l20)/N{0) = 5.3 x 10"^ (Ka/lO' 



-1 ^ - 
cm s 



The atomic oxygen column density along the HD 154368 sight-line is A^(0) = 1.2 x 
10^^ cm~^ (Snow et al. 1996, from absorption line observations of the 1355 A 01] line). 



Based on a search for the C^Bi — X^Ai band of water vapor near 1240 A, ISpaans et al. 
(119981 ) obtained 3 a upper limit on the water column density of 9 x 10^^ cm~^, corresponding 
to A^(H20)/iV(0) = 7.5 X 10~^. This 3 a upper limit lies a factor 70 below the value that 
we would obtain were kra equal to 10~^^ cm^ s~^ as AC06 suggested, and places a 3 a upper 
limit of 1.4 X 10"^''cm^s"^ on kra- 

AC06 suggested that the freeze out of oxygen onto gr ain mantles c o uld d iminish the 
water vapor abundance in molecular clouds, as proposed by iBergin et al.l (120001 ) to explain 
the low H2O abundances measured by SWAS in dense clouds. Ice absorption line observations 
of diffuse/translucent sight-lines, however , indicate that water ice is generally present only 



in clouds of Ay > 3 ( Whittet et al.ll200ll ). Furthermore, in the specific case of HD 154368 



under present consideration, the atomic oxygen is known from direct measurement to be 
1.2 X 10^^ cm"^. The column density of H nuclei along this sight-line, A^h = A^(H) + 2A^(H2 



has been measured to be 4.2 x 10^^ cm~^ ( Snow et al.lll996l ). so the mean line-of-sight oxygen 
abundance is ?7.(0)/nH ~ 3 x 10~^, a value that is ent irely consistent with the abundances 
measured along diffuse sight-lines ( Meyer et al. 19981 ) and inconsistent with a significant 
depletion of oxygen onto ice mantles. 

In summary, the upper limit on the water vapor abundance observed towards HD 154368 
definitively rules out a rate coefficient for the radiative association reaction O+H2 H2O+7 
that is large enough to explain -in the context of the AC06 model- the H2O Iiq— >loi line 
strength measured by SWAS toward IRC+10216. If HSO observations would indicate that 
H2O is formed in the external layers of IRC+10216, an explanation other than the radiative 
association proposed by AC06 would be required to avoid incompatibilities with observations 
toward HD 154368. 



5. Summary 

Our radiative transfer models for H2O in IRC+10216 and other C-rich AGB stars reveal: 

1) H2O in the envelope of IRC+10216 is primarily excited through absorption of photons in 
the 1^2 band at 6 /im and subsequent decay to the ground vibrational state. 

2) The 0-H2O abundance relative to H2 required to account for the 557 GHz lio^ loi 0-H2O 
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line observed in IRC+10216 is in the range (0.5 — 5) x 10 ^, depending on the inner radius 
of the H2O shell. 

3) Multitransition H2O observations with the Herschel Space Observatory will allow us to 
establish the spatial distribution of H2O molecules in lRC+10216, in particular the inner 
radius of the H2O shell, thus discriminating among the different hypothesis proposed for the 
origin (formation or release) of H2O. 

4) A number of other C-rich AGB stars with relatively strong 6 /xm continuum flux are 
expected to be detectable in the lio— >loi 0-H2O hue with HSO, if the H2O outflow rate in 
those sources is similar to that found in IRC+10216. 

5) The relatively low H2O abundance required to explain the lio— >loi 0-H2O emission in 
IRC+10216 makes less restrictive the requirements previously reported for the vaporization 
of icy bodies or Fischer-Tropsch catalysis on the surfaces of metallic grains to work. If H2O 
were formed in the external layers of IRC+10216, we find that an explanation for the H2O 
formation other than the previously reported radiative association O + H2 ^ H2O + 7 is 
required to avoid confiict with observations of the translucent cloud towards HD 154368. 

D.A.N, gratefully acknowledges the support of grant NAG 5 - 13114 from NASA's 
Long Term Space Astrophysics research program. This research has made use of NASA's 
Astrophysics Data System. 
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Table 1. Model parameters for IRC+10216 ^ 

Parameter Value 

Gas velocity field v(r) = 14.5^1 - 0.95^ km s"^ 
Turbulent velocity Vt = 0.65 km s^^ 

H2 density ^(Ha) = c^-^ 

/ 0.72 

Gas temperature Tfe(r) = 12 f ^ j K if ns < 110 

Tk{r) = 10 K if ri5 > 110 
Dust inner radius Ri — 4.5 x 10^^ cm 

Dust outer radius i?„ = 4.0 x 10^^ cm 

Stellar radius i?* = 8.0 x 10^^ cm 

Stellar temperature = 2050 K 

Spectral index /3 = 1.3 if A > 20 /xm 

/3 = 1 if A < 20 /xm 

is the distance to the star in cm, and ri5 = jq^^- 



Table 2. Inner and outer radii for the H2O shell in different models, and derived 0-H2O 

abundances 



Model 


Rint (cm) 


Rout (cm) 


X(o-H20)^ 


A 


4.5 X 10^^ 


4.0 X 10^^ 


4.1 X 10-s 


B 


2.1 X 10^^ 


4.0 X 10^^ 


5.4 X 10-^ 


C 


4.3 X 10^6 


4.0 X 10^^ 


5.1 X 10-^ 



^o-H20 abundances are given relative to H2, 
and an ortho-to-para H2O abundance ratio of 
3:1 is assumed in all models. 
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Table 3. H2O transitions, frequencies, and adopted HSO-HIFI half-power beamwidths and 

sensitivities 

Transition Frequency HSO-HIFI beam'^ Sensitivity'' 

(GHz) (") (10-22 W cm-2) 

lio^loi 556.9 41.0 2.2 

3i2^3o3 1097.4 21.0 12 

3i2^22i 1153.1 18.5 29 

321 ^3i2 1162.9 18.5 30 

523 -^5i4 1410.6 14.5 45 

221 ^2i2 1661.0 14.5 61 

2i2^1oi 1669.9 14.5 62 

3o3^2i2 1716.8 12.5 56 

532^523 1867.7 12.5 65 

2ii^2o2 752.0 29.0 4.1 

2o2^1ii 987.9 21.0 9.2 

lii^Ooo 1113.3 21.0 12 

422 ^4i3 1207.6 18.5 32 

220 ^2n 1228.8 18.5 35 

4i3^4o4 1602.2 14.5 55 

633^624 1762.0 12.5 60 

624-^615 1794.8 12.5 64 



^Adopted from http:/ /www.ipac.caltech.edu/Herschel/hifi/hifi.shtml 



'^Calculated using HSPOT version 1.9.1; sensitivities correspond 
to 5(7 in 1 hr with a spectral resolving power of i? = 10^ 



-24- 




Fig. 1. — Top: Fit to the continuum emission from IRC+10216. Data points are taken from 
Bagnulo et al.l (Il995l ). and the grey hne shows the ISO data. Bottom: Dust temperature 
profile. 
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Fig. 2. — Comparison between the SWAS lio^loi 0-H2O spectrum (from iMelnick et al. 
20011 ) and the hne profile obtained for models A (upper), B (middle), and C (lower), with 
the SWAS beamsize. The spectra of models B and C have been shifted vertically for clarity. 
The required abundances of 0-H2O are given in Table [2l 
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Fig. 3. — Ortho-H20 abundance, relative to H2, required to account for the observed 557 
GHz H2O emission in IRC+10216, as a function of the assumed inner radius of the H2O 
shelL 
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Fig. 4. — a) Fractional populations of the 6 lowest energy levels of ortlio-H20 in model A; 
b) coUisional and radiative pumping rates, per molecule in the Iqi level, for the lio level in 
model A. Solid line: coUisional pumping, which is computed from the Iqi ^ lio coUisional 

rate; dashed line: radiative rate through the path Iqi — > 1^2 lio ~^ lio; dotted line: 

radiative rate through the path Iqi 1^2 2i2 '^21 lio; long-dashed line: radiative rate 
through the path 2i2 —>■ V2 loi — >^ lio- 
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Fig. 5. — Predicted line profile of the 0-H2O lio— >loi line with the HSO-HIFI beam, and 
the telescope axis pointing toward the source center (upper) and toward a position shifted 
41" from the source center (lower). Model A: sohd line; B: dashed hne; C: dotted line. 
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Fig. 6. — Energy level diagram of ortho and para-H20. Arrows indicate the transitions 
whose predicted line fluxes and proflles are shown in Figs. [7] and El 
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Fig. 7. — Predicted H2O line fluxes observable by HSO-HIFI as a function of the inner radius 
of the H2O shell. Adopted beamwidths are listed in Table [3l 
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Fig. 8. — Predicted H2O line profiles observable by HSO-HIFI, for models A (solid), 
(dashed), and C (dotted). Adopted beamwidths are listed in Table [31 
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Fig. 9.— Calculated F(2i2-^loi)/i^(lio^loi) (labelled 1670/557, squares) and 
i^(lii^Ooo)/i^(lio^loi) (labelled 1113/557, triangles) vs. Rm for the IRC+10216 model 
source, but located at a distance of D > 2 kpc (solid lines), D — 0.5 kpc (dashed lines), and 
D = 0.17 kpc (dotted lines). 



- 33 - 



o 
ft 



i 10 



3, 
CO 



-15 



10 



-16 



10" 



10 



-10 



a) Mira: 

Semiregular: 
Irregular: •"^ 



IRC+10216 




b) 




IRC+10216^ 



10" 



10 _^ 10" 

M (Mg yr ^) 



Fig. 10. — a) Continuum flux densities at 6.3 fim, corrected for tlie distance to tlie sources, 
versus tlie mass loss rate for a sample of C-rich AGB stars. Filled symbols indi cate sources 
for w hich the 6.3 /im flux has been directly measured from ISO/SWS spectra (ISloan et al. 
20031 ) . and opened symbols indicate sources for which the 6.3 fim flux has been estimated 
from the available 8.8 fim flux density (see text for details). The grey line shows the fit to 
the general trend achieved with our dust models, b) Same as in a) but divided by the mass 
loss rate. The position of IRC+10216 is indicated in both diagrams. 
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Fig. 11. — Calculated H2O line fluxes vs. M (a), and vs. the 6.3 iim continuum flux (b), 
where the lines are fits of the form F^ne oc -Fg.s/im! results are given for a source at D = 0.5 
kpc. c) Observed continuum fluxes at 6.3 /xm of C-rich AGB stars in function of the distance 
(see Fig. [To] for the meaning of symbols). The lines show the 6.3 /im continuum flux required 
to detect the ground-state lines in 1.5 hours of observing time with HSO-HIFI, if the H2O 
outflow rate is the same as that in IRC+10216 (solid: lio^loi; dotted: ln^Ooo; dashed: 
2i2— ^loi)- Some of the stars that would be detectable in such a case in the lio^loi 0-H2O 
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Fig. 12. — Predicted n(H2)/nH ratio for the best-fit ISpaans et al.l (Il998l ) model, together 
with the ?T,(H20)/?7,(0) ratio that would result if k^a were 10~^^ cm^ s~^. The dashed line 
shows the 3 a upper limit on the observed iV(H20)/A^(0) ratio. 



